Change of an insulator's topological properties by a Hubbard interaction 
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We introduce two dimensional fermionic band models with two orbitals per lattice site, or one 
spinful orbital, and which have a non-zero topological Chern number that can be changed by varying 
the ratio of hopping parameters. A topologically non-trivial insulator is then realized if there is one 
fermion per site. When interactions in the framework of the Hubbard model are introduced, the 
effective hopping parameters are renormalized and the system's topological number can change at 
a certain interaction strength, U — U, smaller than that for the Mott transition. Two different 
situations may then occur: either the anomalous Hall conductivity a xy changes abruptly at U, as 
the system undergoes a transition from one topologically non-trivial insulator to another, or the 
transition is through an anomalous Hall metal, and a xy changes smoothly between two different 
quantized values as U grows. Restoring time-reversal symmetry by adding spin to spinless models, 
the half-filled system becomes a Z2 topological insulator. The topological number v then changes 
at a critical coupling U and the quantized spin Hall response changes abruptly. 

PACS numbers: 71.10Fd, 71.27. +a, 73.43.-f 



Recent interest in non-trivial topological properties of 
insulators [U [5] has spurred intensive search for band 
models with non-trivial topology. This is because of the 
possibility that electron interactions in fractionally filled 
topologically non-trivial bands may lead to the realiza- 
tion of highly correlated fractional quantum Hall states 
[3HB]. On the other hand, cold atomic gases in optical 
lattices with tunable interaction strength open the pos- 
sibility to physically realize topological insulators [8] . 
It is then natural to ask about the effects that electron 
correlations can have on the topological properties for 
proposed models. 

Indeed, non-interacting topological phases are fairly 
well understood, so a great deal of attention has recently 
been given to the effect of interactions [5HT2"] . A local re- 
pulsive interaction, such as that in the Hubbard model, 
drives the system to the Mott insulating state for in- 
creasing interaction strength. Studies of the Kane-Mele- 
Hubbard model [13], H4] have recently been carried out, 
showing that the topological insulator survives until the 
Mott insulating phase is attained [T5RT5] . A topological 
phase can also arise from interactions added to a triv- 
ial band model, leading to a topological Mott insulator 

[T9H22]. 

In the present paper we show that yet another possibil- 
ity exists. Namely, local interactions can drive the sys- 
tem from one topologically non-trivial insulating phase 
into another by changing its topological number, while 
keeping it finite. This mechanism is particularly relevant 
for proposed band models with Chern number C larger 
than one. We present band models which break time- 
reversal symmetry explicitly, a situation analogous to the 
one considered by Haldane for the honeycomb lattice [13] . 
The system is a quantum anomalous Hall insulator when 
the lowest band is filled. If a Hubbard repulsive inter- 



action is also present, the topological Chern number C, 
hence the number of chiral edge states, changes at criti- 
cal values of the Hubbard interaction strength, U, before 
a Mott insulator phase is attained at higher U c > U. 
Such transition is signaled by a change in the quantized 
Hall response. This effect occurs because the interaction 
effectively renormalizes the Hamiltonian parameters for 
the fermions, decreasing the longer ranged hopping with 
respect to the short ranged. Such changes in the effective 
hopping parameters induces a change of the topological 
number of the bands for the fermions. 

The following example models describe electrons in a 
square lattice with two orbitals per site. The Pauli ma- 
trices and (fi — 0,1,2,3) act on the orbital (or 
sub-lattice) space and spin space, respectively, and the 
subscript "0" refers to the identity matrix. The Hamil- 
tonian has the general form: 

H{h) = h(k)-T + h (k)T , (1) 

where h = (h x ,h y ,h z ) and k — (k x ,k y ) denotes the 
momentum vector. The Chern number for the bands 
in Hamiltonian Eq. ([!]) is independent of the choice for 
/io(fe), as computed from the usual expression 

c= lJ dk * dk *w x x w y - h > (2) 

and we therefore neglect ho(k) for the time being, and 
comment on it later. Time-reversal symmetry (TRS) re- 
quires h x ( z ) to be a even function of k and h y to be odd. 
In order to have nonzero C, TRS must not be present. 
The operation of spatial inversion [I : h x (k) — > h x (—k), 
hy( z )(k) — > —h v (z)(—k)] does not change C. 
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FIG. 1: (color online). Spectrum for the model given by 
Eq. |3| in the ribbon geometry: (a) for t% > t[ — 5/4, where 
C = 2; (b) for t 2 < t\ - 5/4, where C = 1. 



i&er. — We consider 



\/2ti (cosk x + cos k y ) , 
\/ 2t± (cos fca; — cos k y ) , 



4t 2 sin fca; sin fc. y + 2i! Y (sin fe^ + sin k y 



k y ) + 5. 



(3) 



The terms t\ and /i y break TRS and are responsible for 
a non-zero Chern number. The terms t 2 and h y break 
spatial inversion symmetry. If t 2 > t[ — 5/4 then Eq. ([2| 
gives C = 2, and if t 2 < t[ - 5/4 then (7=1. The 
spectrum obtained in the ribbon geometry is shown in 
Fig. Qa) for t 2 > t'j - 5/4, the C = 2 case. Two pairs 
of counter propagating edge modes, running along the 
opposite edges of the ribbon, are clearly seen. The case 
t 2 < t'i — 5/4 is shown in Fig. [ljb). Compatible with 
C = 1, now only one edge mode runs along each edge 
of the ribbon. We take the case where the lowest band 
is filled with spin polarized (or spinless) electrons, that 
is, where one electron occupies one of the two orbitals 
per lattice site. The system is then a band insulator 
with a number of chiral edge modes equal to the Chern 
number (see Fig.[T]). We now assume that the parameters 
m Eq. (§ are such that C = 2 and introduce electron- 
electron interactions via the Hubbard term: 



U 



J 3>s 



- 1 



(4) 



where s = 1, 2 denotes the orbital index at lattice site 



J 



and 



"3,8 



is the electron number operator. A suit- 



able approach to the Hubbard model at half-filling is 
the slave-rotor method [23TI25] . In this approach a ro- 
tor e %Bj is assigned to every lattice site and the charge 
at every site is identified with the rotor's angular mo- 



mentum, Li 



n ■, 



1 



-id/dOj. The electron 



operator is decomposed into a fermion and a rotor as 
Cj = fje~ tSj . The interaction can then be simply rewrit- 
ten as: Hi n t — J2j Lj ■ We shall employ the X-boson 
treatment [23] by equating e l6j = Xj with the constraint 
\Xj\ 2 = 1. Such a treatment has the advantage that it 
allows the calculation of the rotor correlation function 
between arbitrary sites of the lattice, as the decay of this 



correlation function with distance in the bose condensed 
phase is of crucial importance in what follows. Previ- 
ous studies with this method concentrated on nearest 
neighbor boson correlation and were focused on the Mott 
transition [TS] . Because the method is well explained 
in Ref. [53] and has been applied to the study of Mott 
transitions in topological insulators a number of times, 
we here simply outline the main steps. The Hamiltonian 
Eq. (JlJ must be Fourier transformed to real space in or- 
der to make the substitution Cj — fje~ l6j and becomes a 
function of /-fermions and X-bosons, H{pX, fX* ). The 
partition function is Z = J D\DfDfDX*DXexp[-S] 
where the action S reads 



l/T 



(It 



X, 



+H(fX, fX*) + ]T »A i (|X,f - 1 )] (5) 
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The imaginary time t is not to be confused with the Pauli 
matrices above. The Hamiltonian in the action reads 



H(fXJX*) = J2T, h i 

ij ss' 



fi,sf 



(6) 



with hij = j^- J2k h(k)e lk '( Ti ~ rj \ where _/V s denotes the 
number of lattice sites. We do the standard mean-field 
decoupling of the fermions and bosons [TH [53] , whereby 
the Hamiltonian becomes a sum of a fermionic and a 
bosonic term, H{fX,fX*) = H f (fJ) + H X {X,X*). 
The /-fermion Hamiltonian, Hf back in momentum 
space is given by expressions in Eq. and Eq. ^ but 
where the hopping parameter between sites i and j is 
now multiplied by {X t X*), i.e., hij -» hu(XiX*). The 
boson Hamiltonian has the form: 



H 



x 



hij(XiX* t 

m: 

Y^h, r ^r"'/,,/,, .v,.\; 

ij ss' 



(7) 



j" 



where j' , j" denote first and second neighbors to the site 
i, respectively, and where A(2) is obtained by averaging 
the fermion fields in Eq. Hx, enjoys in this case the 
square lattice symmetry. The action is accordingly a sum 
of two parts, S — Sf + Sx with 



Sx — 



l/T 



dr 



^ U 

3 



\d T x 3 



3 



i\ 3 (\X 3 \ 2 -l)+H x 



where the rescaling U — > U/2 has been made in order for 
the slave-rotor method to reproduce the correct atomic 
limit [23 , 24 . In going to momentum space and Matsub- 
ara frequency, the boson field can be written as 



x j (r)=^j2 x ^ k y (k ' ri ' VT )+ 



(8) 
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where xq is a density of the condensate of the boson field 
at zero momentum and frequency. The summation in 
Eq. Q excludes the point (iv, k) = (0,0). The action 
for the bosons then reads 



5X 



iv.h 



\X(iv,k)\ 2 



v 

U 



+ Ji7i(fe) + J 2 72(fc) + iA 



(9) 



2cos{k x ) + 2cos(fc y ) and J2(k) — 
Assuming a spatially constant iXj, the 



with 7i(fe) = 
4cos(k x ) cos(fc a ) 
constraint \Xj\ = 1 can be implemented on average in 
imaginary time and space as 



1 = X + 



VU 1 



E 



i 



N s ^ x M7i(fc) + J 2 72(fc)+*A 



(10) 



If the interaction is not strong enough to enter in 
the Mott regime, the condensate density is finite and 
this requires the dispersion relation for the bosons 
yj Ji7i(fe) + J 2 7 2 (fc) + iX to vanish at k = 0. The La- 
grange multiplier is then fixed to the value iX — — 4( Ji + 
J 2 ). The instantaneous boson spatial correlation func- 
tion reads 



(x t xi) = 1 - 



Vu 1 



E 



,ih-(rj-rj 



N s ^ A /J l7l (fe) + J 272 (fe) + iA 



(11) 



The result in Eq. ( 11 ) must be inserted back into the /- 
fermion Hamiltonian Hf, until convergence is attained. 

A simulation with t\ = 1, t 2 = 0.7, t[ = 0.8i 2 and 
5 = 0, yields a transition from the Chern number C = 2 
atJ7<L7toC=lat[/>£7 where U w 1.4. The critical 
interaction for the Mott transition obtained is U c ~ 2.9. 
For U > U we then expect a discontinuous change in the 
Hall conductivity. At U = U the /-fermion bands touch 
and close the gap. As long as t 2 > t\ — 5/4, there is a 
U < U c . If t 2 < t[ - 5/4, then C = 1 for all {/ < f7 c . 

On the Mott insulating side, U > U c , the condensate 
Xo vanishes and the boson correlation function is propor- 
tional to VU, so that the ratio between first and second 
neighbors, (Xi X* +i ) / {XiX* +i+ $) , is independent of U. 
Therefore, the above mechanism is not effective inside 
the Mott phase. 

The physical electron's Green's function has a coherent 
part where the boson condensate xq is the quasi-particle 
weight and the excitations are those of the Hamiltonian 
Hf. We have checked that the Chern number obtained 
from Hf and Eq. ^ agrees with that obtained from the 
physical electron's full Green's function for an interacting 
system (equation (6) in Ref. [H]). In the Green's function 
point of view, the change in the topological number oc- 
curs because a pole of the Green's function moves across 
zero energy [57] . 

Some of the models for flat bands that have been pro- 
posed in the literature, with the purpose of realizing the 



-0.5 



-1.5 



Energy Gap 



C=-2 



AH Metal 
I Hall Conductance 



C=-l 



0.2 




0.6 



FIG. 2: (color online). Topological transition between two 
topologically different regimes through an anomalous Hall 
(metallic) phase as a function of U/U c - The parameters are: 
t = l,ti = 1.5, t% = 1.6, 8 = 0, a = 5. For these parameters 
U c = 12. 6t. The Hall conductance is expressed in units of 
e 2 /h. The energy gap is the indirect gap between the two 
bands. 



fractional anomalous quantum Hall effect, have topolog- 
ical numbers similarly unstable with respect to the Hub- 
bard interaction before the Mott phase is attained. For 
instance, the three-band model presented in Ref. [3] and 
which is expected to be realizable in an optical lattice [8] 
has a transition from C = 1 to C = at U/U c = 0.4, 
according to the above mechanism. 

A model with spin and Z topological number. — We con- 
sider a system with a spinful orbital at each site described 
by a Hamiltonian of the form of Eq. with the matrices 
r replaced by the matrices acting on spin space, a, and 
with 



a sin fc,. 



h„ 



-a sin k x 



hz = 4i 2 cos k x cos k y + 2t\ (cos k x + cos k y ) + 5 (12) 

where h (k) = —2t[cos(k x )+cos(k y )]. Here a is a 
Rashba spin-orbit coupling, the term proportional to t\ 
may be seen as an intrinsic spin-orbit coupling and 5 
as a uniform magnetic field. Time reversal implies in 
this case hi(k) — > —hi(—k). The term h z in Eq. (12) 



breaks TRS. In general, due to the presence of the h a (k) 
term there is an indirect band overlap and the system is 
metallic. Even though the bands are topologically non- 
trivial they will be in general partially filled. Considering 
ho(k) = and placing the chemical potential at zero en- 
ergy one naturally gets a half-filled band insulator with 
non-zero Chern number. The expression in Eq. ^ gives 
C = 2 if I ti| < |t 2 + 5/4\; the Chern number reduces 
to C = 1 if |ti| > |t 2 + 5/4\. Turning on the Hubbard 
interaction, a transition between the different topologi- 
cal phases may be obtained as in the previous model at 
a suitable U < U c . Considering ho, one possible way to 
find an insulating phase is by having a large Rashba term. 
By varying t 2 one finds regimes where the gap becomes 



4 



zero and a band overlap occurs, hence a metallic state. 
Interestingly, increasing further ti takes the system back 
to an insulating phase with different Chern number. The 
width of this metallic phase depends on the Rashba cou- 
pling. Starting from a C = 2 phase and turning on the 
Hubbard interaction we induce a sequence of transitions 
from the topological insulator through an anomalous Hall 
metal and back to another insulating phase with different 
Chern number. This is illustrated in Fig. [2] where we plot 
the indirect gap between the two bands and the Hall con- 
ductance as a function of the Hubbard coupling. There 
is a smooth crossover between the two quatized values of 
the Hall conductance in the two insulating phases. The 
Hall conductance is evaluated using a Kubo formula. 

Model with Z 2 topological number. — We extend the 
model in Eq. ([3]) so as to describe a bilayer square lattice, 
with one orbital per site and where the superscripts in 
t ss denote layer indices. We introduce spin by coupling 
the spin operator o~ z to the terms that break TRS. The 
model reads 



h x = y2ti (cos k x + cos k y ) + t± , 

h y = v2ii (cos k x — cos k y ) o~ z , 

h z — 4i 2 sin h x sin k y + 2t\ (sin k x + sin k y ) o~ z 



(13) 



where t± is real and couples the atom in lattice site j 
of one layer (denoted by an index s) to the one that sits 
directly above it in the other layer s' ^ s. Although 
TRS has been restored by the spin-orbit coupling, the 
model still lacks particle-hole symmetry, however. This 
puts the model in the All class [28] . Studies of the Mott 
insulating phase for systems with spin-orbit coupling and 
TRS, using other techniques, exist in the literature [T5] , 
The half filled system is a Z 2 topological insulator char- 
acterized by a topological number v which is given by the 
parity of the Chern number: v = if C = 2 and v = 1 if 
C = 1. Equation ^ for this case reads 



H(fXJX*) = J2J2 h >- 



fi,s.af 



Y Y* 



where a =t,4- denotes the spin projections. We make 
a similar decoupling of the X-bosons and fermions as 
above, with the difference that the bosons now have an 
extra layer index, Xi S . A calculation with model param- 
eters: ti = 1, t 2 = 0.7, t[ — 0.56, t± = —1 yields a 
transition from v = to v = 1 for U = 0AAU c with 
U c = 4.3. If t\ = 0.42 then one obtains U = 0.84t/ c 
with U c = 4.2. In Fig. [3] we show the phase diagram 
so obtained. Because of TRS, no Hall conductivity ex- 
ists, but the quantized spin Hall conductance, which is 
given by a% y — eC/(2ir), displays an abrupt change at 
the transition, since the number of edge modes changes 
discontinuously (see Fig. [I]). Since v = is not pro- 
tected by additional terms with TRS then, if disorder is 
present, the system may actually go from zero spin-Hall 
conductivity to finite (quantized) spin-Hall conductivity. 
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FIG. 3: (color online). Topological phase diagram for the 
model given in Eq. ( 13 1 with a Hubbard interaction. The 
parameters are: t\ = 1, tz = 0.7, 8 = 0, t± = —1. In the 
phase v — the spin Hall conductance a% y — ej-K whereas 
°xy — e / (27r) in the phase v = 1. 



It is interesting to see that in the Kane-Mele-Hubbard 
model, when either a trivial Dirac mass (staggered poten- 
tial) [13] or a nearest-neighbor Rashba spin orbit coupling 
[29j are present, the present mechanism would destroy 
the topological insulating phase for a certain U < U c . 
This is so because, similarly to the model in Eq 



(131, 



the topological insulator in the Kane-Mele model is in- 
duced by a second nearest neigbor hopping, which is 
thus susceptible to be renormalized by the interaction. 
In graphene, this hopping mimics an intrinsic spin-orbit 
coupling [5^1 [3D] which is already very small. However, 
this mechanism might be relevant for the recently discov- 
ered silicene |31j . where the intrinsic spin-orbit coupling 
is expected to be sizable [52*] . 

Conclusions. — We have shown that a short ranged re- 
pulsive interaction can induce a change from a system's 
non-trivial topological phase to another, and presented 
several examples. The effect occurs because the renor- 
malization of parameters in the Hamiltonian is such that 
the longer ranged hopping is decreased with respect to 
the short ranged. It is an interesting direction for future 
investigations to explore how the renormalization group 
[33) technique for interacting fermions will renormalize 
interaction and kinetic energy parameters so as to induce 
changes in the topological class of a fermionic system. 
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